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Progress in multiloop calculations 

A.G. Grozin‘S 

^'Budker Institute of Nuclear Physics, Novosibirsk 

I briefly summarize the talks on calculation of multiloop Feynman diagrams presented at ACAT’2002 (Moscow 
University). 


Here I briefly summarize the talks about mul¬ 
tiloop calculations iiiiiii which were pre¬ 
sented at the parallel session “Simulations and 
Computations in Theoretical Physics and Phe¬ 
nomenology” at the VIII International Work¬ 
shop on Advanced Computing and Analysis Tech¬ 
niques in Physics Research (ACAT’2002), June 
24-28, 2002, Moscow. 

In order to have some rough classification of 
these contributions, we shall place them on the 
map (Table |^) showing the number of legs and 
loops. The current state of the art is 

• 0 legs — 4 loops (with I mass); 

• 2 legs — 3 loops (massless, HQET, on- 
shell); 

• 3, 4 legs — 2 loops (various cases). 
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1. Kotikov {p2) 

In his talk Q, A.V. Kotikov reviewed several 
methods of multiloop calculations. 


First, he discussed application of the Gegen- 
bauer polynomial technique to the calculation of 
the massless two-loop diagram (Fig. ^ with three 
non-integer powers of denominators. The results 
can be expressed via 3F2 hypergeometric func¬ 
tions of the unit argument with indices tending 
to integers at e —> 0. They can be expanded in 
e in terms of multiple C"Sums (generalizing the 
Riemann function). Reduction of such sums to a 
minimal set of independent ones has been exten¬ 
sively investigated recently. With these results, 
it is not difficult to expand these diagrams up to 
rather high powers of e. 



Figure 1. Massless two-loop self-energy diagram. 


Then he discussed the differential equations 
method. All basis two-loop two-point diagrams 
with a single mass (Fig. where some lines are 
massive) have been calculated using this method. 

2. Grozin {p2, p3) 

In his talk Q, A.G. Grozin reviewed the recent 
progress in calculation of propagator diagrams in 
HQFT and on-shell QCD. 
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It is well-known that there are 3 generic topolo¬ 
gies of massless 3-loop propagator diagrams. 
They can be reduces, using integration by parts, 
to 6 basis integrals. Four basis integrals are triv¬ 
ial. One is the diagram of Fig. with a non¬ 
integer power of the middle line; it is a particu¬ 
lar case of a more general result discussed in the 
previous talk. The last and most difficult basis 
integral (the nonplanar one) can be easily found 
at e = 0 by gluing; there is no easy way to find 
further terms of its e expansion. 

There are 10 generic topologies of 3-loop HQET 
propagator diagrams. They can be reduced, us¬ 
ing integration by parts, to 8 basis integrals. Five 
basis integrals are trivial. Two can be expressed 
via 3 F 2 hypergeometric functions of the unit ar¬ 
gument; they can be easily expanded to a rather 
high power in e. The last and most difficult basis 
integral was found up to the finite term in e, using 
direct integration in the coordinate space. More 
terms of its e-expansion were recently obtained 
using inversion. 

There are 2 generic topologies of 2-loop on-shell 
QCD propagator diagrams with a single non-zero 
mass. They can be reduced, using integration 
by parts, to 3 basis integrals. Two basis inte¬ 
grals are trivial, and the third one is expressed 
via 3 F 2 hypergeometric functions of the unit ar¬ 
gument. However, some of their indices tend to 
half-integers at e —> 0, and the above algorithm of 
expansion in e is not applicable. The case when 
there is another non-zero mass was systematically 
studied recently. There are 4 basis integrals, 2 
of them trivial, and 2 are expressed via 3 F 2 hy¬ 
pergeometric functions of the mass ratio squared. 
Finite parts at e ^ 0 are expressed via diloga¬ 
rithms. 

There are 11 generic topologies of 3-loop on- 
shell propagator diagrams with a single non-zero 
mass. They can be reduced, using integration by 
parts, to 18 basis integrals. The basis integrals 
are mostly known from the calculation of the elec¬ 
tron magnetic moment in QED. Some on-shell di¬ 
agrams are related to HQET ones by inversion of 
Euclidean integration momenta. 


3. Caffo {p2) 

In his talk [^, M. Caffo discussed numerical 
calculation of the two-loop sunrise diagram with 
generic masses (Fig. ||) for arbitrary p^. The fi¬ 
nite parts of the basis integrals satisfy a system 
of linear differential equations. Their values at 
= 0, 00, the threshold (mi -I- m 2 + ma)^, and 
the pseudothresholds like (mi + m 2 — ma)^ are 
known. The differential equations are integrated 
numerically, starting fromp^ = 0, along a contour 
in the lower half-plane. 





Figure 2. The general massive 2-loop sunrise self¬ 
mass diagram. 


4. Kalmykov {p2) 

In his talk Q, M.Yu. Kalmykov discussed pole 
masses of gauge bosons in the Standard Model. 
In order to relate them to the MS masses with 
the 2-loop accuracy, one has to calculate their 1- 
loop self-energies and their derivatives, and the 2- 
loop self-energies, on the mass shell. After expan¬ 
sions in several small parameters, most diagrams 
can be reduced to single-scale on-shell integrals, 
which are expressed via basis integrals using in¬ 
tegration by parts. This is done by the FORM 
package ONSHELL2. For diagrams with light 
fermion loops, the expansion in sin 9\y breaks 
down due to threshold singularities. These di¬ 
agrams are reduced to basis ones by Tarasov’s 
recurrence relations; in their sum, logarithms of 
sin^vv cancel. 
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5. Davydychev {tl, t2) 

In his talk ||], A.I. Davydychev discussed the 
vertex diagrams (Fig. H) on the mass shell = 
M^, = m?). The one-loop diagram can be 

exactly expressed via the one-loop two-point di¬ 
agram with masses M, m and the momentum q 
in 2 — 2e dimensions. The two-loop diagram has 
been calculated, using Mellin-Barnes representa¬ 
tion, at TO = 0 and at non-zero m ^ M. 



Figure 3. One- and two-loop vertex diagrams 


6. Fleischer {tl, bl) 

In his talk [H, J. Fleischer discussed some steps 
towards the two-loop Bhabha scattering calcu¬ 
lation. The matrix element is decomposed into 
Dirac structures times scalar amplitudes. One- 
loop amplitudes, as well as two-loop amplitudes 
with the electron mass counterterm (they are re¬ 
duce to one-loop integrals), have been calculated. 
Tensor integrals are reduced to scalar ones with 
shifted dimension. Scalar integrals are reduced 
to basis ones. One-loop on-shell trianglular dia¬ 
grams are reduced to two-point ones (as in the 
previous talk). 

7. Tentyukov 

In his talk [Q, M. Tentyukov presented the Di¬ 
agram ANAlizer (DIANA). It does not calculate 
diagrams; it does other things necessary for large- 
scale automatic diagram calculations. Diagrams 
for a given process can be generated by QGRAF. 
But calculation of diagrams of each topology must 
follow a separate path. DIANA sorts generated 
diagrams into topologies. Packages for calcula¬ 
tion of Feynman integrals (such as Mincer) re¬ 


quire a specific routing of external and loop mo¬ 
menta for each topology. QGRAF provides some 
momentum routing for each diagram, but not the 
required one. DIANA corrects this deficiency. 
Last but not least, it can produce nice graphical 
representations of diagrams, ready for inclusion 
into a paper. 

REFERENCES 


1 . A.V. Kotikov, hep-ph/0112347 , 

2. A.G. Grozin, these Proceedings, 


ph/0211328 


hep- 


3. M. Gaffo, H. Czyz, E. Remiddi, these Pro¬ 


ceedings, hep-ph/0211171 


4. 


5. 


6 . 


F. Jegerlehner, M.Yu. Kalmykov, O. Veretin, 
Nucl. Phys. B641 (2002) 285. 

A.I. Davydychev, V.A. Smirnov, these Pro¬ 


ceedings, hep-ph/0210171 


J. Fleischer, O.V. Tarasov, T. Riemann 
A. Werthenbach, these Proceedings, hep- 


ph/0210180 


7. M. Tentyukov, J. Fleischer, 
ings, |rep-ph/0210179 . 


these Proceed- 














